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ON THE PRACTICAL APPLICABILITY OF STOKES’ LAW 
OF RESISTANCE, AND THE MODIFICATIONS OF IT 
REQUIRED IN CERTAIN CASES 


By M. S. SMOLUCHOWSKI. 


§ 1. Stokes’ law for the resistance of a sphere in a viscous liquid rests, as is 
well known, on the fundamental assumptions: 


L Slowness of motion, so that the inertia terms in the hydrodynamical 
equations may be neglected, in comparison with the effects of viscosity. 


II. Complete adhesion without slip, of the liquid to the sphere, this being 
considered as a rigid body. 


HI. Unboundedness of the liquid and immobility at infinity. 


In what follows I should like to contribute some remarks on this law with 
regard to certain cases of practical importance, where the underlying conditions are 
changed to some extent, which may be of some interest to those who are ie 
with research work on subjects connected with Stokes’ law. 


First let us touch briefly the question of slipping, connected with the second 
of the above assumptions. Stokes’ calculation can be generalised, by allowing the 
liquid to slip along the surface of the sphere, with a velocity proportional to the 
frictional force in a tangential direction [which in the case of a parallel laminar flow 


3 
implies the surface condition Bu = Hay 


In this case, as Basset has shown, the simple law of Stokes has to be replaced by 


Thus the minimal value of the resistance, for the case of infinite slip (8 = 0), is 
two-thirds of the maximal value for no slip (8 = œ ). 


Now it is generally assumed, on account of the experimental researches of 
Poiseuille, Whetham, Couette, Ladenburg and others, that the slip of liquids along 
solid walls is negligibly small. Mr Arnold’s* recent measurements prove, by their 
exact agreement with Stokes’ law, that the coefficient of sliding friction @ is certainly 
greater than 5,000 and probably greater than 50,000. 


* R. D. Arnold, Phil. Mag. 22, p. 755 (1911). 
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§ 2. On the other side, his experiments, on bubbles of gas moving through 
liquid, gave the unexpected result that the slip at clean* surfaces between gas and 
liquid is infinite, as the velocity turned out too great by 50 per cent. 


Now I think a different explanation of those experiments to be preferable, as 
in the case of gas bubbles or liquid drops also the interior liquid is subject to 
circulation. Some time ago I advised Mr Rybczynski in Lemberg to calculate the 
motion of a viscous sphere through viscous liquid. The calculation is quite easy 
and the result}, published January last year, and deduced also half a year later, quite 
independently of course, by M. Hadamard, is equally simple. It shows that for slow 
motion the inner liquid retains its spherical shape and that the resistance is 

Bp! + Qe 


F= 67 Re un + 3u a E aa a Was (2), 


where w’ designates the viscosity of the liquid in the interior of the sphere. 


Comparison with the above formula shows that the resistance experienced by 
a gas bubble or liquid drop without slip is the same as the resistance of a solid 


, 


sphere with a coefficient of surface friction p=", in fact the velocity and the 


stream lines of the outer liquid are identical in both cases. It would be interesting 
to verify the above formula by experiments on liquids with similar values of u and 
u; in the case of Mr Arnold's experiments the viscosity in the interior was negligible 
in comparison with the viscosity of the outer medium, which had the same effect as 
if the surface slip were infinite. So far his results too are explained without the 
assumption of surface slip. 


§ 3. However, there is a case where the existence of surface slip has been proved 
beyond doubt, namely in rarefied gases. As is well known, the magnitude of the 


coefficient of slipping y =Ë is, aecording to the kinetie theory and also to the old 


8 
experiments of Kundt and Warburg, roughly equal to the mean length of the free 


path of the gas molecules; therefore the phenomenon plays an important part 
even at ordinary pressures in the motion of very minute droplets, as in Millikan’s 
experiments. 


Now unfortunately one cannot use formula (1) for this case, with substitution 
of the empirical value for 8, except for the case of comparatively small slip. For 
if the mean length À is comparable with the dimensions of the moving sphere, the 
ordinary hydrodynamical equations cease to be valid altogether, since the implicit 
assumption underlying them, that the state of the gas is varying little for distances 
comparable with X, is impaired. 


Therefore also the interesting deduction of a corrected formula by Prof. E. 
Cunningham} is not to be considered as a demonstration and Messrs Knudsen and 


* Le, provided the surface be not contaminated with solid films. 
+ W. Rybczynski, Bull. Acad. d. Sciences Cracovie, 1911, p. 40; J. Hadamard, Comptes Rendus, 152, 
p. 1735 (1911); 154, p. 109 (1912). 
t E. Cunningham, Proc. Roy. Soc. 83, p. 357 (1910). 
M. G. II, 12 
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S. Weber may be right in trying to get closer approximation by other, purely 
empirical formulas*. At any rate the formula proposed by Cunningham 


-1 
W o eiie Ë Ae 
R 
serves remarkably well for interpolation, considering the experiments of those authors 
and those of Mr McKeehanf. It is preferable to write it in the form 


B 
‘lies 
w 
where p is the density of the gas, as mistakes are easily involved by using the mean 
length of free path A, which is a very indefinite term and really has no precise 


meaning. 


F=6rpRe 


> 


For great rarefaction the resistance is proportional to the cross-section of the 
sphere, and for this case the calculation can be carried out exactly, if the way is 
known, how the interaction between the surface of the sphere and the gas molecules 
takes place. If they rebound like elastic bodies, we get in accordance with 


Cunningham 
4 ae : 
F = 3 a R Tr pc Ve 


where V is the square root of the mean square of molecular velocity. The numerical 
coefficient, as calculated from the experiments mentioned above, is considerably 
larger, it amounts to 1°65 (Knudsen and Weber) or 1°84 (McKeehan). McKeehan 
concludes that molecules are reflected from the surface of the sphere only in a normal 
direction; I think, however, that his theoretical formula is not quite exact and at 
any rate his conclusion seems to me at variance with fundamental principles of the 
kinetic theory of gases. I think that the experimental results are explained best by 
the view, supported also by other researches of this kind, especially those of Knudsen, 
that a solid surface acts in scattering the impinging molecules irregularly in all 
directions whether with or without change of mean kinetic energy. We shall not 
go into these questions now, however, as they belong to the kinetic theory of gases, 
not to hydrodynamics. 


§ 4. Now let us consider what modifications are required in Stokes’ law, if the 
third of the fundamental assumptions is impaired, the liquid being limited by solid 
walls, or a greater number of similar spherical bodies being contained in it. 


In this case the linear form of the hydrodynamical equations makes it possible 
to attain their solution by a method of successive approximations, analogous to the 
method of images used in the theory of electrostatic potential. It consists in the 
successive superposition of solutions formed as if the fluid would extend to infinity, 
but so chosen as to annul the residual motion at the boundaries, with increasing 
approximation. 

This method was used first by H. Lorentz in order to determine the influence 
of an infinite plane wall on the progressive movement of a sphere, and we shall refer 


* M. Knudsen and 8. Weber, Ann. d. Phys. 36, p. 981 (1911). 
{| McKeehan, Physik, Zeitsch, 12, p. 707 (1911). 
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to his formulas later on*. He found that the resistance of the sphere is increased 
° 


: : 9 š Š , 
by a fraction amounting to + for normal motion for parallel motion, if @ 


9 R 
> 16a 
denotes the distance from the wall. Mr Stock in Lemberg has extended the calcula- 
tion for the second case to the fourth order of approximation, including terms 


: JENS 
with (=) T 
In a somewhat similar way Ladenburg} calculated the resistance experienced by 


a sphere, when moving along the axis of an unlimited cylindrical tube, and his result, 
indicating an increase in comparison with the usual formula of Stokes in the propor- 


tion of 1: 14247 (where p=radius of the tube), has been verified with very 
satisfactory approximation by his own experiments and by those of Mr Arnold. 

§ 5. Now let us apply this method to the case where a greater number of 
similar spheres are in motion, and extend a little further now an investigation which 
I had begun in a paper published last year§$. Imagine a sphere of radius R, moving 


with the velocity c along the X-axis, its centre being situated at the distance z from 
the origin. It would produce at the point P (with coordinates #, n, £) certain current 


velocities Uo, Vo, W, of order =, defined by Stokes’ equations, if the fluid be unlimited. 


But if we assume this point P to be the centre of a solid sphere of radius R, we have 
to superpose a fluid motion u,v,w,, chosen so as to annul the velocities of the primary 
motion at the points of this sphere and satisfying the conditions of rest for infinity. 


This motion may be called the “reflected” motion; it can be found with any 
degree of approximation, by making use of the solution of the hydrodynamical 
equations given by Lamb, in form of a development in spherical harmonics. But 


as it is of order w: at the surface of the second sphere, which is its origin, it 
seems probable, a priori, that its magnitude at the first sphere will be of order 


c (y, and I have verified this as well as the following results by explicit calculation. 


Thus if we confine ourselves to terms of order (P), we can apply a simplified 


method of evaluating the mutual influence of such spheres, by neglecting the 
difference between the velocity at the centre of the second sphere and at its surface. 


That is to say, the sphere P, being at rest, is subjected to frictional forces 


X =6rpRu, 
Y=6rphy, 
4=6rphu, 


* H. A. Lorentz, Abhandlungen ü. th. Physik, 1. p. 23 (1906). In Millikan’s determinations of the ionic 
charge the increase of resistance due to the presence of the condenser plates may produce an increase of 
the order of one-thousandth. 

+ J. Stock, Bull. Acad. d. Sciences Cracovie, 1911, p. 18. 

t R. Ladenburg, Ann. d. Phys. 23, p. 447 (1907). 

$ M. 8. Smoluchowski, Bull. Acad. d. Sciences Cracovie, 1911, p. 28. 
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on account of the motion of the first sphere; on the other side, the moving sphere 
experiences a reaction by virtue of the presence of the sphere P, such as if this 
would execute simultaneously the three motions — to, — %, — w; the three current 
systems resulting therefrom, according to the usual formulas of Stokes, produce at 
the centre of the first sphere nine current components, giving rise to nine com- 
ponents of frictional force, to be calculated each according to Stokes’ law of resistance. 


If both spheres are in simultaneous motion, the mechanical effects are found by 
superposition of the forces corresponding to the two cases where one of them is 
moving and the other one at rest. 


In this way an interesting conclusion is obtained for the case where both 
spheres are moving in parallel directions with equal velocity: then both are sub- 
jected to equal additional forces in the same direction, one component in the direc- 


tion of motion, tending to diminish the resistance by the amount ete |: _ Ee "| 
the other component along the line joining the centres, towards the sphere which 
is going ahead, of amount = hı —_ iç] [where @ is the angle between 


the line of centres and the direetion of motion]. 


Thus two heavy spheres of this kind would sink faster than Stokes’ law is 
indicating and, besides, their path must be deflected from the vertical towards the 
line of centres by an angle e defined by 
| 3 R 


sin e= —— Wami 


es |sin 8 cos 8, 


§ 6. Analogous methods are applicable to a greater assemblage of spheres. 
The motion results from superposition of simpler solutions, where one sphere is 
supposed moving and all the other ones resting. Each of the component solutions 
comprises the direct action, and for higher approximation also its “ reflections.” 


Now if the parallel motion of a cloud of n similar spheres is considered, the 
resistance of each of them will be diminished by an expression proceeding after 
powers of R, the first term of which will be of the order of magnitude poRP t, We 
see that these developments would be divergent for an infinite number of spheres. 
It is evident that for instance an infinite row of spherical particles, arranged at equal 
distances, would acquire infinite velocity, by virtue of their gravity, as also an infinite 
cylinder would behave in the same way. This applies a fortiort to two-dimensional 
infinite assemblages. 


Stokes’ law of resistance will not be true even approximately, and the develop- 


i : nR. 
ment will cease to be convergent in general, unless —— is small, where S denotes 


S 


a kind of mean distance, comparable with the linear dimensions of the cloud. 


§7. The same result follows from the following simple reasoning. Imagine 
a spherical cloud of radius SŠ, containing n spherical particles, each of radius R and 
density =, suspended in a medium of viscosity p, of negligible density, for example 
a cloud of minute drops of water in air. Then currents will take place in the 
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spherical cloud and it will attain a certain velocity as a whole, which may be 
calculated after the formula (2), just as if the cloud would form a homogencous 


3 
medium of density n (4) and of the same viscosity as the outer medium. The 
L 


4 w ga 


mass velocity resulting therefrom, of amount „s — , is superposed upon the 


15 Su 
displacement of the particles, relative to the moving cloud, taking place with 
velocity s Thus evidently the downward velocity will be much inereased, and 


nR 


Stokes’ law cannot be true even approximately, unless 5 is small in comparison to 


unity. This condition shows that Stokes’ law can be applied only to particles con- 
stituting clouds of exceedingly scarce crowding, and it is easily seen that it would be 
quite erroneous to apply it to actual fogs or actual clouds in the atmosphere, with 
diminished transparency [as in this case the aggregate cross-section of the particles 
nR? is comparable with the cross-section of the cloud Sr]. As an illustration how 


: : ' 0 Be 
cautious we must be in this respect, I may mention that the ratio ae amounts to 


8 
M” 
10 and even to 100 for a cubic centimetre cloud as produced by Sir J. J. Thomson 
and H. A. Wilson, in their experiments on the determination of the ionic charge. 


§ 8. What has been said applies of course only to clouds moving in an other- 
wise unlimited medium. The conditions of motion are quite different for a cloud 
contained in a closed vessel, as in the experiments just referred to. Prof. E. 
Cunningham has attempted to evaluate the order of magnitude of the correction 
to be applied to Stokes’ law in this case. His estimate is founded on the supposition 
that each particle moves approximately in such a way, as if it were contained in a 
rigid spherical envelope, of radius comparable with half the distance to its next 
neighbours. Now this supposition does not seem quite evident, although we shall 
see that it leads to a result of the right order. 

We can calculate the resultant motion in quite an exact way, if we consider 
a homogeneous assemblage of equal spherical particles, moving all of them with the 
same velocity c in the direction of negative X, towards an infinite rigid wall, which 
we assume to be the plane YZ. In this case we see, by making use of H. A. Lorentz’s 
calculation before alluded to, that a moving sphere x,y,z produces at a point £, 
situated on the axis of X, a velocity component 

aa 2Fela E52) ] 4 BM A Ole 
; 4p P p 

The first part of this expression, containing r= V(ze—Ë) ++ 2, is the com- 
ponent of direct motion, according to Stokes; the second part is the component 
caused by “reflection” at the plane YZ; it contains the distance between the point 
Ë and the reflected source p= V(z + EP + 3⁄2 + 2. 


A e e R i 
The terms with higher powers of = have been neglected, as we confine ourselves 


to the first approximation. The total current produced in the point £ by the motion 
of all the particles is equal to U = Xu, where the summation is to be extended over 
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all their values of x, y,2 Now we might consider it right to replace the summation 
by an integration, since one particle corresponds to a space A’, if A denotes a sort 
of mean distance between the particles. In this case the result would be very 


simple, for we should have 
1 
U= ga] [| vdedyae, 


The integrals of the separate terms constituting « can be evaluated explicitly 
if we extend them to a cylinder with YZ as basis, of height h and of radius Z. Then 
we can use the well-known expression for the potential of a disk in points of its axis, 
and expressions derivable from it by differentiation with respect to & and by these 
means we find the unexpected result that the integral current U is zero, if we extend 
the summation to an infinite value of G. 


But in reality U is not defined by integration but by summation. Evidently 
both operations lead to the same result for distant parts of the space, but not for 
those parts whose distance from the point £ is comparable with the distances A 
between two particles. Therefore the resultant current U in points at a great 
distance (in comparison with A) from the wall will be given by 


3 Re 
U=+ 8 } 


l 

eren (4), 
l as A añ | 

where B= +]J + A da du de — sf (1 + A J 


to be extended over a space great in comparison with Æ, is a purely numerical 
coefficient. 

In order to evaluate 8 we must know how the particles are arranged. If we 
suppose an arrangement in rectangular order, we can get easily an approximate value 
by explicit calculation and by integrating over a cube of height H, constructed 
around the point & which gives 


(fll a ; 
{I (a + S) dedyda = 8H" Hog +v3)—$log2— 7]. 


A 
2 
sion for 8 is rapidly converging with extension of the limits of integration. In this 
way I have found the approximate value @= 3:09, and therefore the resistance for 
one particle will be 


It is sufficient to take H equal to a small uneven multiple of “5, as the expres- 


1+ — = baphe 


Paes | ancien 
TA 1+ 23 ; j (5) 
This formula would apply, of course, also if the particles were arranged in a different 
way, but then the numerical value of 8 would be different. Our result agrees to the 
order of magnitude with Prof. Cunningham’s estimate, which led him for the case 


F = 6rpRe 


! me š 3 R. 
of an equilateral arrangement to a similar formula, with a coefficient of 1 included 


within the limits 3°67 and 45. 


§ 9. However, the practical application of this formula is rather questionable, 
as it applies only to a regular arrangement of particles. If they were arranged in 
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clusters, the correction might even become negative. It is interesting to note that 
the average value of £, for a particle whose position relatively to the other ones is 
defined by pure accident, would be zero, and that seems quite natural, as the average 
current of liquid U in the cross-section must be zero. Thus it follows, what we 
should not have expected at first sight, that Stokes’ law applies for the particles of an 
actual cloud, on an average with no correction whatever, of this order of magnitude. 


The evaluation of the quadratic terms would be much more complicated of 
course, as then all possible kinds of single reflections caused by any one sphere have 
to be taken into account. 


The general result of our calculation shows at any rate that Stokes’ law is 
undergoing but small corrections if applied to the particles of a uniform cloud filling 
a closed vessel. But it is important to note that things will change entirely if the 
cloud is not of quite uniform density or if it does not fill the whole empty space 
between the walls. Then as a rule convective currents will arise, which in certain 
cases may be of preponderant influence. Their velocity may be calculated ap- 
proximately by considering the medium as a homogeneous liquid subjected to 
certain forces, the intensity of which per unit volume corresponds to the aggregate 
force acting on the particles contained in it. 


Consider for instance an electrolyte in an electric field. If it is conducting in 
accordance with Ohm’s law, the average electric density is zero and no currents will 
take place. But in bad liquid conductors, with deviations from Boyle’s law, con- 
vective currents may arise, which may influence also materially the apparent value 
of the conductivity. They have been observed long ago, for instance by Warburg*. 


Similar movements may be produced in ionised gases, and I think more atten- 
tion ought to be paid to them than usually is done. In experiments where the 
saturation current of strong radio-active material is observed between condenser 
plates wide apart+, these phenomena may be of importance as producing an ap- , 
parently greater mobility of the ions than under normal conditions. 


§ 10. There is another application of the theoretical methods exposed above 
which may be mentioned. Imagine a two-dimensional infinite assemblage of equal 
spherical particles, distributed uniformly over the plane æ =}, whilst the plane YZ 
again may be supposed to be a rigid wall. Now let all these particles be moving 
along the plane in direction Y with equal velocity c; what motion will be produced in 
the surrounding liquid, and what will be the resistance experienced by every particle ? 


According to Lorentz again the motion produced by a single sphere moving 


parallel to a fixed wall is, when higher powers of the ratio = which we suppose 


to be a small quantity, are neglected : 
= PRE £ “j 2 . oxy? (ar 
pad ke 1+(2) _ 3 Re i4 (5) a a aa 
4r x] 4p p/ | 2 p° p 


where the first term is the direct current according to Stokes, while the remaining 
terms represent the current reflected by the wall, just as in the former example. 


bo! 


* E. Warburg, Wied. Ann. d. Phys. 54, p. 396 (1895). 
T Cf. Rutherford, Radio-activity, pp. 35, 84. 
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We might also in this case calculate the resultant current by forming Sv over 
all values of y and z, and derive therefrom the resistance of a single particle. But 
we shall confine ourselves to the following remarks. 


In the extreme case where the particles are so crowded, as nearly to touch one 
another, a lamellar flow will take place in the liquid between the fixed wall and the 
plane æ=} with a velocity =°, while on the other side of the plane æ=} the 


liquid will be dragged along by the sheet of moving particles with the constant 
velocity c. The frictional force per unit of surface of the plane æ=} is evidently 


equal to Ë © therefore the resistance experienced by each particle is 


l 


which is much smaller than Stokes’ law would indicate, as A is of the order of R but 
the distance J is supposed to be of higher order. 


Now consider the other extreme case, where the distances A between the 
particles are so great that Stokes’ law is approximately valid, which requires A to 
be of order l. Let us calculate the resultant motion of the liquid for points at 
infinite distance from the wall (ë = œ). For such points the summation mentioned 

à 1 x 1) ie ee} il 
above can be replaced by integration; besides we can put ee: F = “3 Z= 2 = 6i ; 
m ; a 
and thus we get 
; 9 Relé f ydi 
Ë: = S= a || ues > 
(& re 7° + 22) 


This integral can be transformed by putting y= s sin $, z= s cus ó$, dydz = sdsdó, 
and we get finally 


6 Rize 
Teas 
x A! 
By comparing this with Stokes’ law for the resistance # we have 
V. = =a ; ; 
Ag 


that means that in both cases the liquid at a great distance from the wall will be 
dragged along, in a parallel direction to it, with such a velocity as if the force 


corresponding to unit surface were distributed uniformly over the liquid, in a 


F 
A? 
plane at a distance / from the fixed wall. This result, which can be gencralised for 
a greater number of similar layers, seems natural enough if the distances between 
the particles are small in comparison with their distance from the wall, so that the 
assemblage can be considered as if forming a homogeneous medium, but we see it 
remains true for particles widely apart. Without going into further details, I may 
only mention that this result has an important bearing on the theory of electric 
endosmose, which will be explained elsewhere with full details. 


§ 11. I may conclude with a brief remark about the influence of the inertia 
terms in the hydrodynamical equations (assumption I), which have been neglected 
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as well in Stokes’ original calculations as in the above reasonings. It is well known 
that this neglection is justified only if the ratio Sor is small in comparison to unity. 
H 


But it has been proved by Oseen* in an important paper, commented upon in a very 
interesting way by H. Lamb, that the solution given by Stokes is defective even 


if this criterion is fulfilled; for at distances r where 7“ is large, the inertia terms 
H 


must be of prevalent influence over viscosity. Oseen himself has given a solution 
which is different from Stokes’ equations for those distant parts of the space and 
gives better approximation there. However, the resistance of the sphere depends 
only on the state of movement in its immediate neighbourhood, therefore the resist- 
ance law of Stokes is not impaired by those results. The condition of its validity 
may be defined more exactly by means of the recent experiments of Mr Arnold, 
which have shown that it holds with very good accuracy (one half per cent.) for 
spheres moving under influence of gravity, provided their radius is smaller than 


0:6 z, where the critical radius 7 is defined by the relation T =]1. This means that 


? ` 


the ratio fog must be smaller than (0'6¥ = 0:22. 


§ 12. The inertia terms are of greater importance, in the case before alluded to, 
where the motion of a greater number of similar spheres is considered. For it is 
legitimate to calculate the forces of reaction between such spheres by using Stokes’ 

` equations for slow motion only if they are lying within the space where viscosity 
is predominant over inertia. Mr Oseen has generalised recently the calculation 
of the interaction of two spheres given by me by introducing in it his solution of 
Stokes’ problem. The forces exerted on the two spheres come out unequal in this 
case and are given by much more complicated expressions. They become identical 
with the first approximation given by me if the distance r between the two spheres 


satisfies the condition that i is small. Mr Oseen thinks this to be a great restric- 
be 


tion on the validity of those formulas for experimental purposes, but he omits the 
factor o in the above expression. We satisfy ourselves easily that, for instance, in 
the case of water-drops in air, as in Sir J. J. Thomson’s and H. A. Wilson’s con- 
densation experiments, the limit of validity for r is of the order of several centi- 
metres; in Perrin’s experiments on the applicability of Stokes’ law to the particles 
of emulsions it would amount to hundreds of metres. It is also sufficiently great 
for direct experiments, when highly viscous liquids are used, as Ladenburg did in his 
elaborate research. Ordinary hydraulic experiments, with water and spheres of a 
size to be handled conveniently, are excluded of course when Stokes’ law or any of 
those modifications are in question. 


One might try to apply Oseen’s method of approximate correction for inertia 
also to the other cases treated above, but it will imply rather cumbersome calculations 
and, besides, for movements in closed vessels it will be generally of lesser importance 
than in a liquid extending to infinity. 

* Oseen, Arkiv f. mat. astr. fysik, 6 (1911); H. Lamb, Phil. Mag. 21, p. 112 (1911). 
+ F. Oseen, Arkiv f. mat. astr. fysik, 7 (1912). 


